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ABSTRACT 

We use the Sloan Digital Sky Survey, along with the NRAO VLA Sky Survey and 
the Westerbork Northern Sky Survey to define a sample of 746 radio-loud quasars and 
measure their 330 MHz to 1.4 GHz spectral indexes. Following previous authors we 
take the radio spectral index as an indicator of the orientation towards the quasars 
such that more pole-on sources tend to have flatter spectral indexes. We use this proxy 
for the orientation of quasars to investigate the effect orientation may have on optical 
spectra. 

Quasars with flatter spectral indices tend to be brighter. However, we find no in- 
dication of reddening in steep-spectrum QSOs to indicate obscuration of the accretion 
disk by a torus as a possible explanation. Nor do we find increased redddening in the 
flat-spectrum sources which could imply a contribution from jet-related synchrotron 
emission. 

We reproduce a previously-described anti-correlation between the width of the 
Mgn line and radio spectral index that indicates a disk- like geometry for the Mgn 
BLR. However, in contrast to previous authors we find no such correlation for the C iv 
line suggesting a more isotropic high-ionisation BLR. 

Both the [O n] and [O in] narrow lines have more flux in steep spectrum sources 
while the [O in] /[On] flux ratio is lower in these sources. To describe both of these 
effects we propose a simple geometric model in which the NLR exists primarily on the 
surface of optically thick clouds facing the active nucleus and the NLR is stratified 
such that higher-ionisation lines are found preferentially closer to the nucleus. 

Quantitatively we find that orientation may effect the observed strength of narrow 
lines, as well as ratios between lines, by a factor of ^2. These findings have implications 
for the use of [O in] and [O n] emission lines to estimate bolometric luminosities, as 
well as comparisons between narrow line luminosity functions for type 1 and type 2 
objects and the potential of emission-line diagnostic diagrams as an accurate tool with 
which to distinguish types of active galactic nuclei. 

Finally we find no evidence that BAL QSOs have a different spectral index dis- 
tribution to non-BALS although we only have 25 obvious BALs in our sample. 
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1 INTRODUCTION 

Orientation plays a key role in the unificatio n model for ac- 
tive galactic nuclei (AGN; lAntonuccll Il993) as one of the 
most important parameters that define the optical charac- 
teristics of an AGN. While we have become familiar with 
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appealing to orientation in order to explain the different 
classes of observed AGN, the effects of orientation within a 
single class of objects are less well understood. This paper 
aims to study the effect of orientation on spectra of quasi- 
stellar objects (QSOs) by using the radio spectral index as 
a proxy for orientation. Our findings can broadly be cate- 
gorised into those pertaining to the broad-line region (BLR) 
and those pertaining to the narrow-line region (NLR). We 
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will discuss these separately throughout much of the paper. 
We also touch on the effect of orientation on the optical 
continuum shape and proportion of broad-absorption line 
(BAL) systems in QSO spectra later in the paper. 



1.1 Orientation and the BLR 

The BLR is too small to be spatially resolved with cur- 
rent telescopes, hence the geometry of the BLR, and the 
physics that governs its dynamics, is not well constrained. 
There are some indications that the dynamics of the BLR 
may be dominated by the local gravitational potential (e.g. 
iPeterson fc Wandell Il999h . Virial supermassive black-hole 
(SMBH) mass estimators make use of this assumption to 
estimate the masses of SMBHs in QSOs. 

Integral to SMBH mass estimators for QSOs is an as- 
sumption about the geometry of the BLR. This is commonly 
expressed as a single factor, /, which translates the mea- 
sured FWHM of a spectral line into the velocity of the BLR, 

«BLR, 



«blr = v7 x FWHM. 



(1) 



Hence, given the radius of the BLR (Yblr), and assuming 
that the local gravitational field governs the dynamics, the 
SMBH mass (Mblr) can be estimated as 



Mb 



/ 



tblrFWHM 2 
G 



(2) 



The value of / is unknown. Model values give / = 3/4 for 
a spherically symmetric BLR l|Peterson fc W andcl 1993), or 
/ = 1/(4 s in 2 8) for a disk BLR incl ined at an angle 8 to the 
observer (|McLure fc Dunlod 12001 ), hence orientation can 
affect the value of / a nd can bias SMBH mass estimates 
I Jarvis fc McLurdl2002h . 

Given the interest in BLR dynamics, and their impor- 
tance to SMBH mass estimation, there have been many 
studies of the structure of the BLR. Techniques employed in- 
clude defining a proxy fo r orientation and using it to look at 
the e ffect on broad lines dBrotherto n 1996; J arvis fc McLurd 
2006), working out SMBH masses through a method inde- 
pendent of the virial technique and then working backwards 
through the virial equation to give f IIMcLure fc Dunlod 
120021 ; lOnken et all |2004| ; lLabita et al-lboOrJ ). looking at the 
broad line distribution as a whole and glean ing results 
from that (|Osterbrockl 1 19771 ; iFine et al l IMPOST), velocity- 
resolved reverberation mapping dKollatschnv fc Dietrich! 
1 19961 ; iDennev. et al.|[201ol : 
through X-ray absorption 



Bentz. et al.l 20101) and recently 



Risaliti et aill2010l ). 



In radio-loud QSOs the ratio (R) between core and lobe 
emission has been employed as an orientation indicator un- 
der the assumption that core emission is D oppler-boosted 
when observed down the ax i s of the object (|Orr fc Brownd 
1 19821 ; I Wills fc Browed 1 19861 ; IVestergaard et al.ll2000l) . Core 
emission from radio-loud QSOs tends to have a flatter spec- 
tral index (a > —0.5; S v oc v a ) than lobe emission due to 
a superposition of many Doppler-boosted synchrotron self- 
absorbed spectra. Hence observing a radio loud QSO down 
the jet axis, not only enhances the o bserved core emission , 
but also flattens the radio sp ectrum (I Wills fc Browndl 19861 ; 
I Jarvis fc McLurdl2002l 120061 ) . 

The relation between a and observing angle is by no 




e 

Figure 1. The 330 MHz to 1.4 GHz radio spectral index as a 
function of source orientation from the empirical simulation of 
Wilman et. al. (2008). 



means exact and must include significant scatter due to in- 
trinsic differences between radio sources. A guide for con- 
necting a to the observing an gle for a radio-loud Q SO can 



be found in the simulations of Wilman, et al. (2008) based 



on the dual-population models of Jackson fc Walll (1 19991 ). In 
thei r semi-empirical simul ation of the radio source popula- 
tion IWilman. et alj (|2008l ) included a relativistic beaming 
model for sources with a random distributio n of viewing an- 
gles o n the sky. Full details can be found in lWilman. et all 
(2008), and in Fig. □ we plot the 330 MHz to 1.4 GHz radio 
spectral index against the viewing angle to the source for all 
FR I and FR II objects in the simulation. Note that in these 
models the extended emission is assumed to have a = 0.75 
and the resulting distribution of spectral indexes is limited 
by the same value. The solid line in Fig. [T] is a linear fit to 
the data in log(a + 0.75)-# space. We only plot the figure for 
8 < 45° as objects observed from a larger angle would likely 
be obscured. Steeper-spectrum objects do exist and help il- 
lustrate the difficulties in modelling the radio source pop- 
ulation. However, the radio spectral index has been shown 
to work as an or ientation indicator (I Wills fc Brownd 1 19861 : 
iBrothertonll 19961 ). though we must be mindful that there is 
a large scatter in the relation. 

By taking R or a as proxies for orientation, and com- 
paring these with broad line profiles from optical spec- 
tra, several authors have found strong evidence that broad 
emission lines are narrower in radio-loud QSOs viewed 
down the jet axis dWills fc Browndll986l; lBrothertonlll996l ; 
IVestergaard et al.ll2000l ; I Jarvis fc McLurdl2006l ). The simple 
interpretation of these results in that the BLR is confined to 
a disk, hence reducing the line-of-sight velocity when viewed 
down the pole. It is also worth pointing out that if the BLR 
was in any way coupled to the jet, i.e. the material was en- 
trained by the jet, then we would expect a large offset and 
broadening of the emission lines which we do not observe. 

For radio-quiet QSOs, a common technique for study- 
ing the effect of orientation on the BLR is by making 
orientation-independent empirical estimates for the mass of 
QSO SMBHs. By comparing these with virial SMBH mass 
estimates, which would be effected by a non-spherical BLR, 
/ can be calculated and used to deduce the emitting struc- 
ture. 

iMcLure fc Dunlod (|2002h compared SMBH mass esti- 
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mates from the M — a relation and the H/3 virial relation 
in a small sample of radio-quiet QSOs. Based on a fit to 
their line width distribution, whe re they assume a pu r e-disk 
model for the BLR (i.e. «» =0), iMcLure fc Dunlod (|2002h 
calculated the average correction factor / for virial SMBH 
mass estimates as a function of observed line width. They 
then showed that correcting virial masses by this factor im- 
proved the agreement with SMBH mass estimates based on 
the M — a relation. Hence they argue that their disk model 
for the H/3 BLR is valid. 

Since, lLabita et all (|2006h and iDecarli et all (|2008l ) 
have produced similar analyses indicating that the C IV BLR 
may also be disk-like. 

Studies of the broad line distribution on it s own have 
also s hed light on the geometry of the BLR. lOsterbrockl 
(1977) argued that the lack of narrow-line Seyfert Is in a 
sample of 36 objects was evidence that rotation c ould not be 
invoke d to explain the width of broad lines. Later lFine. et all 
(2008) and lFine et"al] (|2010h investigated the Mgn and C iv 
line width distribution in a large dataset of (~ 100, 000) 
QSO spectra. These analyses showed that the variation in 
line width between QSOs is so slight that the BLR cannot 
be s trongly geometr i cally c onstrained. 

iRichards. et~afl (2002) showed that C IV emission-line 
shifts with respect to other spectral features are due 
to/accompanied by a drop off in flux in the red wing of the 
line. They proposed an orientation-based model to describe 
their observations in which the C IV BLR is associated with 
an outflow, the far side of which can be obscured depending 
on the angle to the source. 

A somewhat stable empirical picture of the BLR is 
emerging from these observational studies. Each of the most 
often studied broad lines (H/3, Mgn and Civ) shows indica- 
tions that the basic geometry may be disk-like. However, the 
Mg n and C IV emission lines at least (and likely H/3) cannot 
be confined completely to a disk and must have some other 
component to their velocity field. 

1.2 Orientation and the NLR. 

The NLR can be observed in type 2 AGN, and many stud- 
ies into the effect of orientation on NLR properties exploit 
this by comparing type 1 and type 2 objects. The implica- 
tion being that type 1 AGN are all viewed roughly pole-on, 
and type 2 AGN are viewed roughly edge-on. Therefore any 
differences between the NLR properties in suitably matched 
samples of type Is and type 2s can be attributed to orien- 
tation. We note here that this paper is concerned primarily 
with high luminosity QSOs and the role of luminosity on 
the torus in AGN and the critical angle between type 1 and 
type 2 is not f ully known (although see lWillott et al]|2000t 
ISimpsonl feOQS). Hence results on low-luminosity Seyferts 
may not directly translate to quasars. 

Studies comparing the flux ratio between high and low- 
excitation lines indicate that t ype 1 AGN have a higher- 
ionisation NLR than type 2s Isiiuder fc Osterbroc k 198 it 
ICohenll 19831 ; ISchmitd[l99Sl ; iMuravama fc Taniguchilll998bl ). 
A number of interpretations for this result have been sug- 
gested. The most popular calls for a stratified NLR in which 
higher-excitation lines are located in a smaller region and 
can be hidden by the obscuring tor us when viewed edge on 
jMuravama fc TaniguchHll998bl lah. ISchmittl (|l998l ) also ar- 



gued for a second interpretation in which type 1 AGN have 
smaller NLRs, and hence are ionised by a stronger contin- 
uum. 

Subsequent studies, eit her comparing the differences be- 



tween type 1 and 2 spectra (Nagao et al.il2 00l| ; Zhang et al 



2008), or spatially resolvin g the NLR (jBennert et al 
l2006bl lal; iKraemer et al.l 120091 ) tend to support a stratified 
NLR in which higher-excitation lines are obscured when 
viewed edge on. However, little work has been done in as- 
sessing the impact of orientation in large samples of type 1 
quasars. 



1.3 This paper 

In this paper we fol low on from work presented in 
Ijarvis fc McLurej (|2006l ) and use the radio spectral index 
as a proxy for orientation to investigate the effects of orien- 
tation on both the broad- and narrow-line region in radio- 
loud quasars. In Section [2] we present the sample used in 
this analysis, Section [3] describes the emission line analysis 
that is new to this work and in Section!?] we present our re- 
sults and discuss them in terms of models for the broad and 
narrow-line regions. Throughout, this paper we assume a flat 
(n m ,n A ) = (0.3,0.7), H = 70kms" 1 Mpc" 1 cosmology. 



2 THE SAMPLE 

Both the photometric and sp ectroscopic optica l data for this 
paper co me from the SPSS (I York, et all 2000) data release 



5 CDR5; lAdelman-McCarthv et al 



compiled by ISchneider. et all (|2007l ). Details of the SPSS 



20071 ) quasar catalogue 



tele scope and spectro g raph are given in lGunnT et all (|2006l ) 
and lStoughton. et all |2002l ). 

In creating the SPSS QSO sample used in this pa- 
per, ISchneider. et~ail (|2007f ) visually inspect all of the can- 
didate spectra to determine their classification. Hence the 
final sample includes all broad-line SPSS spectra, combin- 
ing radio and optically selected objects in amongst other 
broad-line objects selected under different criteria. To en- 
sure we are not biased, we include in our sample only ob- 
jects which were colour-selected a s low or high-redshift QSO 
candidates (|Richards. et al.ll2002T ) as the BEST target. And 
to ensure accurate radio identification we take only those 
objects which lie in the F aint Images of the Radio Sky 
(FIRST; iBecker et all fl 9951 ) region, and have a FIRST de- 
tection (6255 objects). 

This core sampl e was then cros s mat ched with the 
WENSS catalogue (|Rengelink et all Il997l ) to produce 
746 matches. Thes e were then matched to the NVSS 
|Condon et al.l fl998) for which we found all sources had a 
unique match. 

The result is a sample of broad-line galaxies from the 
SPSS with WENSS and NVSS fluxes for each. Note that, 
while all of the WENSS detections have an NVSS detection 
there are many NVSS sources in the cross-matching area 
with no WENSS detection, due to the relative depths of the 
surveys. Furthermore the requirement of a broad line in the 
spectrum means that contamination from BL Lac objects, 
that may be strongly beamed in the optical, is low. In fact 
taking the classic broad-line equivalent width cut of 5 A (e.g. 
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Figure 2. (a) The spectral index distribution of our sample as calculated from the 330 MHz and 1.4 GHz fluxes. (6) The relation between 
radio spectral index and r-band absolute magnitude. 



lUrrv fc Padovani| [l995) to define BL Lacs there are none in 
our sample. 

Fig. [2]» shows the 330 MHz to 1.4 GHz radio spectral 
index distribution of our sample. In this paper we are taking 
radio spectral index as an orientation indicator. The objects 
with flatter spectral indices are considered to be viewed more 
pole-on than the source s with steeper radio sp ectra. Fig. &{b) 
shows how the r-band l|Fukugita et al.|[l996l ) absolute mag- 
nitude of our sample depends on radio spectral index. The 
correlation (a Spearman rank test gives r s ~ —0.23 with 
P{fa) <C 1%) evident between a and M r may indicate a 
degree of obscuration of the accretion disk in QSOs viewed 
from an angle away from the pole (however see Fig. [TU] and 
discussion later in paper). Furthermore the correlation in 
Fig. [2f 6) may be indicative of a selection bias in our sample. 
Since we find brighter objects have steeper spectral indexes 
the optical selection of the SDSS sample may bias us towards 
these (potentially more pole-on) sources. In the analysis that 
follows we will try to correct for potential luminosity bias. 



3 EMISSION LINE ANALYSIS 

In this paper we will compare a series of emission line param- 
eters with radio s pectral index to inf er con clusions ab o ut the 
emitting regions. iFine. et all (120081 ) and iFine et al.l (|2010r i 
have already performed an analysis of all of the Mgn and 
C IV line in our sample and we will use the parameters for 
these lines derived in those papers. In addition we will be 
measuring parameters of the narrow [O n] A3727 and [O ill] 
A5007 lines. 

[O n] is not blended with any other emission lines and 
is therefore relatively straightforward to analyse. We fit a 
quadratic continuum between windows either side of the line 
and subtract it from the spectrum. The resolution of the 
spectra does not require us to fit the [O n] doublet as two 
individual emission lines, therefore a single Gaussian is fitted 
to the line to define its extent. We then sum the flux in the 
continuum-subtracted spectrum within ±1.5 X the FWHM 
of the fitted Gaussian. 

[O ill] is a doublet, and can be blended with the wing of 
the broad H/3 line. We fit a continuum between the extrem- 
ities of the H/3 line and subtract it from the spectrum. We 
then fit a double- Gaussian model to H/3 and single Gaus- 
sians to the [Om] A4960 and 5007 doublet fixing their rel- 
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Figure 3. Example of a fit for the [O III] line for QSO J000149.52- 
100251.7. The top plot shows the original spectrum and linear 
continuum fit. The second plot shows the continuum subtracted 
spectrum along with Gaussians fitted to the H/3 and [O in] A4960 
lines (dashed), and the [Om] A5007 line (solid). Vertical lines 
indicate the ±1.5FWHM limits within which we sum the flux. 
The y-axis in each plot gives the flux density normalised to the 
peak in the top plot. 



ative amplitude to 1:3. The three Gaussians not pertaining 
to [O ill] A5007 are subtracted from the data and the flux 
was summed between ±1.5x the FWHM of the Gaussian fit 
to the line. 

Note that this method for fitting the [O ill] lines does 
not always give accurate fits to the whole H/3 region of the 
spectrum. Differing H/3 profiles, continuum shapes and vary- 
ing iron emission in the region all have an effect. However, 
out fitting method does accurately correct for contaminants 
to the [O ill] A5007 line. However, to make sure the results 
were acceptable, all fits were inspected manually and in some 
cases the continuum fit had to be adjusted to obtain ade- 
quate fits. An example fit is shown in Fig. [3] 



4 RESULTS 

Here we compare various emission line parameters with the 
radio properties of our sample in order to derive conclusions 
about the line emitting regions in QSOs. Since the BLR 
and NLR differ in size by many orders of magnitude we will 
consider them separately. 



4.1 The BLR 

We are interested in how broad lines in quasar spectra are 
effected by orientation. In particular, previous studies have 
shown evidence that emission lines are broader in objects 
observed edge on, indicating a disk-like BLR. We are using 
the radio (330 MHz to 1.4 GHz) spectral index as a proxy for 
orientation in our sample, and in Fig. [4] we show Mgn (a) 
and Civ (6) emission line widths plotted against the radio 
spectral index. 

In each panel of Fig. U there is a significant amount of 
scatter. However, it is apparent that there is an anticorre- 
lation between Mgn line width and a in Fig. |4ja), most 
evident for a < —0.2. A Spearman rank test on the data in 

Fig. Eta) gives r s 0.26 (P(\r 3 \ > 0.26) < 0.01% given 

a random distribution). On the other hand Fig. |4j6) shows 
no evidence for a correlation (r s ~ —0.02; P(r 3 ) ~ 0.79). 
There are fewer objects with a C IV line width measurement 
and we test whether this is reducing the significance of any 
correlation with respect to Mgn. We resample the Mgn line 
width distribution randomly to have the same number of 
objects as have C IV measurements and perform Spearman 
rank tests in each of the subsamples. In all cases we found 
P(r s ) < 0.05 and P(r s ) < 0.01 in 80% of the subsamples. 

To highlight how the profile of the lines depends on 
a, Fig.[5]shows median- value (variance- weighted composites 
show the same results) composite spectra of the Mgn (a) 
and Civ (b) line regions for objects that have a < —0.5 
(solid) and a > —0.5 (dashed). In each figure the spectra 
have had a continuum subtracted between the red and blue 
ends of the plot, and have had their flux densities normalised 
to their maximum values for comparison. 

The correlation in Fig. |4ja) reproduces the results of 
(|Jarvis fc McLure! 120061 ) and is in agreement with other 
studies indicating a flattened BLR for the Mgn line. The 
lack of any correlation in Fig. |4j6) may be of more conse- 
quence. The indication is that the C IV BLR is, to a large ex- 
tent, isotropic and unaffected by or ientation effects. This is 
in contrast to prev ious studies (e.g. IVestergaard et al]|2000l ; 
iDecarli et a"l]|2008h . 



4-1.1 Discussion 

Our results for the Mgn line are in line with previous re- 
sults. The implication is that the Mgn emission region has 
a significant velocity component in the plane of the disk, po- 
tentially orbital motions. However, while the difference be- 
tween the composite spectra in Fig. [5] is highly significant, 
the magnitude of the difference is not huge. This indicates 
that a thin-disk model for the BLR is not sufficient and some 
other, non-planar, component to the motion is required. 

We find no evidence for orientation effects in the 
C IV line. Previous authors have reported an effect (e.g. 



Orientation effects in quasar spectra 5 




O I ' 1 ' 1 ' 1 ' 1 

2000 4000 6000 8000 

FWHM (km/s) 

Figure 6. Here we compare the disk model from McLure &; Dun- 
lop (2002), that connects the observed FWHM of a line to the av- 
erage correction factor / for a virial SMBH mass estimate (solid 
line; note McLure &; Dunlop (2002) define / as the square root 
of the parameter in equation [2] we plot \ff to present a simpler 
comparison to their data), to a model in which the FWHM of 
lines is unrelated to orientation or SMBH mass (dashed line). The 
two models show almost identical behaviour indicating that data 
which follow the trend of these lines is not necessarily evidence 
for a disk-like BLR. 

IVestergaard et al ] |2000l : IDecarli et aljliooih and we discuss 
the potential reasons f o r the discrepancy here. 

IVestergaard et all (|2000l ) compared the radio core-to- 
lobe flux ratio R with C IV line widths for 37 radio-loud 
QSOs. They found no correlation between the FWHM of the 
C IV line and R. However, they found that line width mea- 
surements that were more weighted towards the wings of the 
line (the full width at 20 % maximum, and inter-percentile 
velocity, IPV, widths) did show a signific ant anti-correlation 
with R. Given the small sample size in IVestergaard et all 
(2000), and the difficulties they had in finding a correlation, 
it may be that th ei r resu lt is not robust. 

IDecarli et alj (120081) follow a method laid out by 
iMcLure fc Dunlod (|2002r ) in which virial SMBH mass es- 
timates are compar ed to empirical mass est imates from the 
Magorrian relation (Ma gorrian. et aD[l9 98). The argument 
follows that residuals between the two mass estimates can 
be ascrib ed to orientation eff e cts. A ssuming a disk model for 
the BLR IMcLure fc Dunlod (|2Q02T] derive the average cor- 
rection f actor for virial mass estim ates (/ from equation [2] 
note that lMcLure fc Dun lop 2002 define / as the square root 
of / as we define it in equation [5} as a function of observed 
FWHM. Put simply, in their model an emission line that is 
observed to be narrow is more likely to be observed pole-on, 
and hence effected strongly by orientation, than an emission 
line that is observed to be broad. Hence narrower lines will, 
in general, have a higher value of ^/f (= 1/(2 sin 8) in their 
model). 

Both IMcLure fc Dunlod (|2002h and IDecarli et al] 
( 2008) find that this disk model for the BLR give good fits to 
the observed residuals between virial SMBH mass estimates 
and those from galaxy scaling relations. Furthermore, the 
disk model gives an improved fit compared to an isotropic 
model in which / = const. However, we argue here that this 
is not necessarily evidence for a purely disk-like BLR. 

In recent years some concern has grown up over the 
accuracy of viri al SMBH mass estimators. In particular 
iFine et all (|2010h argued that there is not enough dynamic 
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Figure 4. These figures show Mgn (a) and Civ (b) line widths as a function of radio spectral index, a. While there is considerable 
scatter in (a) there is a trend such that flatter spectrum sources have narrower line widths. (6) shows no evidence for a trend (see text 
for details). 
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Figure 5. Composite spectra of the Mgn (a) and Civ (6) line regions for objects with a < —0.5 (solid) and a > —0.5 (dashed). In (a) 
there is a clear difference between the two lines with the flat spectrum sources giving a narrower emission line. In (b) the Civ line shows 
very little difference between the flat and steep spectrum sources. 



range in either Mg II or C IV line widths for them to be use- 
ful in virial estimators given the typical uncertainty of the 
mass estimates. If we assume that emission line width does 
not help in mass estimation, and is essentially random, then 
narrower lines will bias mass estimates low, and broad lines 
will bias mass estimates high. The result, in terms of derived 
values for / as a function of observed FWHM, is extrem ely 
similar to the disk model of iMcLure fe Dunlod d2002l). In 
Fig. [5 ] we show the best fit model from IMcLure fc Dunlod 
(2002) which assumes V or b = 4375 km/s, a or b = 1400 km/s 
and 6 max = 47° as the solid line. We have then taken 
a model in which observed line widths are unrelated to 
the virial mass. T h at is, we have taken the same Vorb as 
IMcLure fc Dunlod (120021 ) and assumed that the FWHM 
is independent of the virial velocity of the BLR, hence, 
J(FWHM) = Vorb/FWHM. This relation is plotted in 
Fig. [6] as the dashed line. Note that our line is not a fit to 
data, we simply take parameters from the lMcLure fc Dunlod 
(2002) fit which assumes a different model. However, the two 
lines show the same tren d. Hence, we argue that the resid- 
uals iDecarli et alJ (120081 ) find between their SMBH mass 
estimates may be due more to a failing in the virial mass 
estimators than a disk-like BLR. 

There may be some concern that the Mg 11 and C IV 
lines occur in spectra of QSOs in differing redshift ranges, 
and hence different luminosity ranges. However, studies of 



the widths of these emission lines in large samples of QSOs 
have shown little-to-no dependence on luminosity or red shift 
jCorbett, et al.ll2003l ; lFine, et alJl20oj ; lFine et al.ll2010r i. As 
a further test we restrict our analysis to only those (88) 
objects in our sample that have both Mg 11 and C IV in their 
spectra. Correlating the line widths of these with their radio 
spectral index we still find a significant correlation (r s ~ 
-0.3; P(r 8 ) < 0.01) for Mgn and none for (r s ~ -0.1; 
P{r s ) ~ 0.2) Civ. 



4.2 The NLR 

The ratio between the flux of the [O ill] A5007 and [O n] 
A3727 lines gives a measure of the ionisation level of the 
NLR. One of the primary drivers of this is the intrinsic lu- 
minosity of the source and in Fig. [7] we show how the ratio 
depends on the r-band luminosity of the QSO. 

Previous studies of the NLR have found that high exci- 
tation lines are emitted from a smaller region when com- 
pared to low excitation lines. Furthermore, when viewed 
edge-on high-excitation lines are less luminous, potentially 
due to absorption close to the nucleus. To test the effect 
of orientation on the NLR in our sample of QSOs we plot 
the [O in] to [0 11] flux ratio against radio spectral index in 
Fig. Ufa). To check that the results are not due to luminos- 
ity bias we correct for the correlation evident in Fig. [7] by 
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Figure 7. The flux ratio between the [O in] and [O n] lines plotted 
against the r-band absolute magnitude of the object. The clear 
correlation indicates that the NLR in brighter sources is in a more 
highly ionised state. The solid line is a minimum x 2 fit- 



Figure 10. g — r colour of our sample correlated with radio spec- 
tral index. There is no evidence for a correlation here, indicating 
that QSOs viewed from an angle are not heavily reddened com- 
pared to those viewed pole-on. 



dividing the measured flux ratio by the value of the best 
fit relation at that magnitude. These corrected results are 
shown in Fig. [8^6). 

At face value the correlation apparent in Fig. [8] is ex- 
actly as we expect. For fiat spectrum pole-on QSOs there 
is less absorption of the inner high-excitation NLR and the 
[O ill] /[On] ratio is increased. However, if we look at the 
fluxes of [O n] and [O ill] individually we find that the picture 
is more complicated. Fig.[9]shows the (a)[Oll] and (6)[Om] 
line luminosities plotted against the spectral index of the 
QSO. 

In both Fig. [9] (a) and (6) the same behaviour is ob- 
served. There is an anticorrelation between the line flux and 
spectral index implying that, on average, more pole-on QSOs 
have fainter oxygen lines. Taking a = —0.5 as the approxi- 
mate mid-point of the spectral index distribution, those ob- 
jects with steep-spectra (a < —0.5) have ~1.3 times the 
[O in] flux and ~1.4 times the [O n] flux compared with the 
flat-spectrum (a > —0.5) objects. This is not easy to square 
with our interpretation of Fig. [8] above: that the reduced 
oxygen flux ratio for steep spectrum sources is due to ab- 
sorption. 

One possibility is that the flat-spectrum QSOs are in- 
trinsically fainter in the optical, and therefore produce fewer 
ionising photons. While Fig. [2] shows that the observed op- 
tical magnitudes are, in fact, brighter for the flat spectrum 
QSOs. It may be that in the steep-spectrum QSOs, that are 
observed from a larger angle, may suffer from more extinc- 
tion in the optical. Extinction should also redden the optical 
spectrum of our sample, and in Fig [10] we show the g — r 
colour of our sample as a function of radio spectral index. 

We find no evidence for a correlation in Fig.[l0](a Spear- 
man rank test gives r s ~ —0.03 with P{r a ) ~ 0.4). Since the 
optical bands sample different rest wavelengths at different 
redshift we check that differing redshift distributions for low 
and high a were not responsible for the lack of correlation 
by also measuring the optical spectral index off the SDSS 
spectra. This is an inherent part of our line fitting proce- 
dure and we find that, regardless of which line we measure 
the optical spectral index under, we find no correlation with 
radio spectral index. It is also worth noting that we also 
find no evidence for increased reddening in the flat-spectrum 




Figure 11. A schematic showing the NLR model discussed in the 
text. The central source ionises the face of optically thick clouds 
in the surrounding galaxy. In the inner NLR the radiation field is 
stronger and higher ionisation lines are brighter compared to the 
outer NLR. When observed from an angle a large portion of the 
ionised surface of each cloud is visible and the lines are bright, 
however, the inner high-excitation NLR is obscured by the torus. 
When observed from the pole the inner NLR is not hidden by the 
torus, but we are looking through the back of the NLR clouds 
and see less line emission in total. 

population, which could arise due to synchrotron emission 
associated with the jet. 

The implication is that the QSOs in our sample that are 
viewed from an angle are not heavily reddened compared to 
those viewed pole-on. Hence we find it difficult to explain 
Figs. [8] and [9] through dust extinction. 

Another potential explanation for Figs. [5] and [5] is with 
a model in which the narrow-line emitting clouds of the NLR 
are themselves optically thick. In this model only the surface 
of the clouds facing the central ionising source would be 
ionised and, when observed from behind, the ionised region 
would not be visible through the cloud. Fig. [11] shows a 
schematic of this model. 

In the model presented in Fig. llll the NLR exists on the 
face of optically thick clouds ionised by the central source. 
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Figure 8. The flux ratio between the [Om] and [On] lines plotted against the radio spectral index for the QSOs in our sample. There 
is a clear positive correlation indicating that the NLRs of more pole-on QSOs have a higher excitation level. There is an indication of a 
bimodal distribution in a in this plot with an increased number of objects with a ~ 0. This is likely due to Doppler boosting of sources 
from the highly-populated fainter end of the luminosity function. In (b) we correct for potential luminosity bias using the best fit line in 
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Figure 9. The luminosity of the [Oil] A3727 (a) and [Om] A5007 (f>) lines plotted against the spectral index of the source. For spectra 
in which the emission line was not visible an upper limit is plotted calculated from the spectral S/N. Both (a) and (6) show the same 
behaviour with line flux anti-correlated with a. This is primarily due to a lack of objects with flat spectral indices, i.e. large a, and bright 
emission lines. 



In the inner NLR, the radiation field is stronger, and higher 
ionisation lines are brighter compared to the outer NLR. 
When observed from an angle a large portion of the ionised 
surface of each cloud is visible and the lines are bright, how- 
ever, the inner high-excitation NLR is obscured by the torus. 
When observed from the pole the inner NLR is not hid- 
den by the torus, but we are looking through the back of 
the NLR clouds and see less line emission in total. This 
model can simultaneously explain the relation between the 
[O ill] /[On] ratio and radio spectral index, and the trend 
between narrow-line flux and a. 



4-2.1 Discussion 

In the previous section we have shown that the optical 
narrow-line strengths of QSOs correlate with their radio 
spectral index that we employ as an orientation indica- 
tor. Narrow line strengths are commo nly employed to esti- 
mate bolometr ic luminosities of AGN (|Heckman et al.ll2004l ; 
ISimpson|[2005h . as well as being exten sively used in diagnos - 
tic analyses such as BPT diagrams (|Baldwin et al.lll98ll ). 
Here we discuss the implications of our findings with respect 
to these analyses. 



To quantify the potential effects of orientation we split 
our sample by the rough midpoint of the radio spectral in- 
dex distribution (a = —0.5). We find that the differences in 
average line luminosity between the low and high a bins are 
~ 0.17 dex and ~ 0.12 dex for [On] and [Oin] respectively. 
These values give the approximate lower bound for the effect 
of orientation on observed line strengths. 

In fact, the effect of orientation is likely to be consid- 
erably more than quoted above since the considerable scat- 
ter in the relationship between a and QSO orientation (e.g. 
Fig- HI nas the effect of smoothing the observed relationship. 
Quantifying the bias introduced through taking a as a proxy 
for orientation is difficult due to the many unknowns. How- 
ever, we can estimate the magnitude of the effect through a 
simple simulation. 

We start by producing a rectangular radio spectral in- 
dex distribution with equal probability for — 1 < a ns < 0.. 
We use the subscript ns to indicate that for this original a 
distribution we are assuming there is no scatter in the a— 
orientation relationship. We then produce an observed radio 
spectral index (a bs) distribution by simply adding random 
Gaussian noise to these values at the level found in Fig. [T] 
(rms~ 0.25). We also produce narrow- line luminosities from 
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the a„„s with a linear relation, including Gaussian scatter, 
in which the gradient, zero-point and level of scatter are 
varied so that the resulting a (, s -line luminosity correlation 
resembles those observed in Fig(9] 

Splitting our simulated sample by a bs = —0.5 and 
a ns = —0.5 and comparing the average line luminosity in 
the low and high a bins for each split we can obtain an esti- 
mate for the level of bias in the results quoted above. In the 
case of [Oil], where we found a 0.17dex offset between the 
bins separated by a j, s = —0.5, we estimate the unbiased 
offset to be ~ 0.25 dex. In the case of [O III], where we found 
a 0.12 dex offset between the bins separated by a b s = —0.5, 
we estimate the unbiased offset to be ~ 0.20 dex. That is, 
[On] lines are roughly twice as strong (and [Om] lines ~1.6 
times as strong) in QSOs observed from an angle, when com- 
pared with those observed pole on. 

[O in] line luminosity has been employed to estimate 
the bolometric luminosity of AGN, m ost often in type-2 ob- 
jects in which the co re is obscured (jHeckman et al.l |2004| ; 
lLamastra et al.l 120091 '). However, the bolometric corrections 
must be calibrated against unobscured type-1 objects. Ex- 
trapolating our results from QSOs to type-2s, it may be that 
orientation effects result in the [O in] luminosities of type-2 
objects being > 1.6 times the value of objects observed polc- 
on. The result would be that bolometric luminosities esti- 
mated for these objects would be similarly o verestimated, as 
well a s Ed dington ratios in the ex amples of iHeckman et all 
(|2004 ) and lLamastra et all l|200gh . 

Our results also impact comparisons between [O ill] lu- 
minosity funct i ons between t ype-1 and type-2 AGN (e.g. 
ISimpsonl 120051 ; iReves, et al.l |200S|). and the c onclusions 
drawn from these comparisons. " ISimpsonl |2005l ) find the 
type-1 [Om] luminosity function is consistently lower than 
the type-2 luminosity function in their sample of SDSS 
AGN. However, if we shift their type-1 luminosity function 
to higher luminosities by a factor of 1.6 we find that, for 
[Om] luminosities > 2 — 3 x 10 34 W the type-1 luminosity 
function is higher than the type-2. This then impacts their 
estimation of the type-l/type-2 AGN fraction, increasing it 
significantly. 

Since we find that the ratio between [O n] and [O m] 
line luminosities is also affected by orientation, it may be 
that the accuracy of BPT diagrams is limited by orientation 
effects. We find the [Oii]/[Oiii] flux ratio rises by at least 
a factor of two with increasing a. [O n] has a comparable 
ionisation potential to Hn, and so we may expect similar 
variability in the [O m]/H/3 ratio: the y-axis on the standard 
BPT plot. Hence more pole-on sources would be biased low 
in standard BPT plots, making them appear more like star- 
forming galaxies. 

4.3 BALs 

The nature of broad absorption systems in QSOs is still 
under debate. A common model assumes that all QSOs 
have a BAL outflow, most likely accelerated equatorially, 
and only in those obj ects observed through the outflow 
do we see a BAL (e.g. IWevmann et al.lll99lT ). Spectropo- 
lametric analysis have provided evidence to back up this 
picture l|Cohen et al.lll995l ; iGoodrich fc MiMerlll995l L How- 
ever, more recent studies using radio properties to gauge 
the viewing angle to the source have complicated the issue. 




a 



Figure 12. The radio spectral index distribution of Civ BAL 
QSOs and C IV non-BALs. The height of the shaded histogram 
has been multiplied by two to highlight the shape of the BAL dis- 
tribution. We find little difference between the two distributions 

I Becker et al ] l|200d ) showed that radio selected BAL quasars 
can display compact radio morphologies and possess both 
steep and flat r adio spectra indicating no preferen tial orien- 
tation (see also lZhou et~ai]|2006l ; lboi' et alj|2009h . 

There are relatively few Mg H BALs, while C IV exhibits 
BALs in ~ 20 % of objects (|Trump. et al.ll2006T ). We there- 
fore look back at all the C IV lines in our sample and cat- 
egorise them as BAL or non-BAL objects by eye. Of 240 
QSOs with the C IV line in the spectrum we find 25 objects 
which display definite BALs and 171 which definitely do not. 
In the remaining 44 objects there is some indication of ab- 
sorption but the S/N is such that we cannot be sure of a 
BAL. 

Fig. [12] shows the spectral index distributions of the 
definite BAL and non-BAL samples, to highlight the shape 
of the BAL distribution the height of the shaded histogram 
has been multiplied by two. We find that these distribu- 
tions are not significantly different with respect to a KS test 
(Dks = 0.17; P(Dks) = 0.5), and our results are not altered 
if we include the objects which have indefinite BAL identifi- 
cations. However, our numbers are small here and we cannot 
definitely say whether the orientation orientation distribu- 
tion of BAL vs non-BAL objects is the same. 



5 CONCLUSIONS 

We use the SDSS, NVSS and WENSS to define a sample 
of radio-loud QSOs and measure their 330 MHz to 1.4 GHz 
spectral indexes. Using the radio spectral index as a proxy 
for orientation towards the QSOs in our sample we investi- 
gate the effect of orientation on their optical spectra. 

In our sample the Mgn line is significantly broader for 
objects with steep spectral indexes indicating disk-like ve- 
locities for the Mg II BLR. However, we find no such correla- 
tion for the C IV line potentially indicating a more isotropic 
high-ionisation BLR. 

Both the [O n] and [O ill] narrow lines have more flux 
in steep spectrum sources while the [O ill] /[On] flux ratio is 
lower in these sources. To describe both of these effects we 
propose a simple geometric model in which the NLR exists 
primarily on the surface of optically thick clouds facing the 
active nucleus and the NLR is stratified such that higher- 
ionisation lines are found preferentially closer to the nucleus. 
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We discuss the implications of our finding on analyses that 
employ narrow line fluxes (e.g. estimating bolometric lumi- 
nosities/BPT diagrams) and show that orientation can skew 
results by as much as a factor of two. 

We find that QSOs with flatter spectral indexes tend 
to be brighter. However, we find no indication of redden- 
ing in steep-spectrum QSOs to indicate dust obscuration of 
the accretion disk as a possible explanation, or in the flat- 
spectrum population which would imply contamination by 
synchrotron emission from the jet. Finally we find no evi- 
dence that BAL QSOs have a different spectral index dis- 
tribution to non-BALS although we only have 25 obvious 
BALs in our sample. 
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